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ABSTRACT 



Exoplanet transit and Doppler surveys discover many binary stars during their op- 
eration that can be used to conduct a variety of ancillary science. Specifically, eclipsing 
binary stars can be used to study the stellar mass-radius relationship and to test predic- 
tions of theoretical stellar evolution models. By cross-referencing 24 binary stars found 
in the MARVELS Pilot Project with SuperWASP photometry, we find two new eclips- 
ing binaries, TYC 0272-00458-1 and TYC 1422-01328-1, which we use as case studies to 
develop a general approach to eclipsing binaries in survey data. TYC 0272-00458-1 is 
a single-lined spectroscopic binary for which we calculate a mass of the secondary and 
radii for both components using reasonable constraints on the primary mass through 
several different techniques. For a primary mass of Mi = 0.92 it 0.1 Mq, we find 
M2 = 0.610 ± 0.036 Mq, Ri = 0.932 ± 0.076 Rq and R2 = 0.559 ± 0.102 Rq, and find 
that both stars have masses and radii consistent with model predictions. TYC 1422- 
01328-1 is a triple-component system for which we can directly measure the masses and 
radii of the eclipsing pair. We find that the eclipsing pair consists of an evolved primary 
star (Ml = 1.163 ± 0.034 Mq, Ri = 2.063 ± 0.058 Rq) and a G-type dwarf secondary 
(M2 = 0.905 ±0.067 Mq, R2 = 0.887± 0.037 Rq). We provide the framework necessary 
to apply this analysis to much larger datasets. 



1. Introduction 

Photometric and Doppler radial velocity (RV) surveys have generated comprehensive data sets 
in their quest for extrasolar planets. These surveys have been successful in discovering several hun- 
dred planet^ and furthering our understanding of how planets form, their dynamically evolution, 
and their distribution of physical properties. These data sets also provide a wealth of ancillary 
science studies involving variable stars, binary stellar systems and brown dwarf companions. The 
amplitudes of the photometric or Doppler variability from such objects are often several orders of 
magnitude larger than extrasolar planet signals, making them comparatively simpler to identify 
and characterize. 

Cross-referencing stars in both photometric and Doppler surv eys can minimize t he amount of 



follow-up data required to conduct these ancillary projects (e.g., iKane et al. II2009I ). For binary 
systems, there are several advantages in combining data from both types of surveys. One disadvan- 
tage of using spectroscopic binaries (SBs) as a source of detached eclipsing binary (EB) candidates 
is that the line-of-sight inclination of the orbital plane is unconstrained. The a priori probability 
that a given SB eclipses can be estimated as {Ri + R2) /a, where Ri and R2 are the radii of the 
primary and secondary, respectively, and a is the orbital semimajor axis. Typical a priori values 



^http : //exoplanet . eu/| 
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range from a few percent up to ~ 25%, therefore conducting photometric fohow-up results in null 
detections a majority of the time. 

In addition, the predicted ephemerides for eclipses based on modeling the Doppler RV data 
depend sensitively on the orbital phase coverage, the precise values of the eccentricity and argument 
of periastron, and the amount of time elapsed since the last measurement. Conversely, archival 
data from photometric surveys consist of tens of thousands of epochs over several years, resulting 
in nearly continuous phase coverage for short orbital periods. Given the orbital period from the 
spectroscopic orbit, it is more efficient to check archival photometry from transit surveys than to 
attempt follow-up photometric observations based on the transit ephemeris from the radial velocity 
data. 

Attacking the problem from the other direction results in a different set of drawbacks. The 
advantage of starting with an EB detected from a photometric survey is that the favorable geometry 
for an eclipse is already known to exist. However, obtaining high-precision, time-series RV follow- 
up is a resource-intensive enterprise. Such observations are conducted using heavily-subscribed 
instruments (generally Echelle spectrographs) on moderately large telescopes (usually larger than 
two meters in diameter) for which access to observing time is quite competitive. Obtaining the 
best orbital solution nominally requires a minimum of seven epochs that sample the orbital phase 
well, requiring many nights of follow-up observations. 

Doppler RV surveys, on the other hand, already have the necessary number of epochs and 
sampling rate because an identical observing strategy is required to detect planets. The observing 
cadences are selected with great care in order to provide synoptic coverage of all possible exoplanet 
orbital periods and phases. Some surveys cease observations if a star exhibits Doppler variation of 
many km over short timescales, indicati ve of an SB, to m aximize planet yield. Other surveys. 



such as the multi-object MARVELS survey (jGe et al. II2009I ). continue to observe the SB to fully 



characterize the orbit. For the specific case of MARVELS, the multiplicity of a given target is usually 
not known prior to the start of observations. Continuing to observe SBs to fully characterize their 
orbit works particularly well in multi-object surveys, where keeping a few SBs per field maximizes 
science results while not significantly decreasing the capability of detecting extrasolar planets. 

Another disadvantage of EBs that only have photometric survey data is that the masses of the 
two components can only be estimate d. Although the rnass of the primary can be estimated using 



colors (e.g., the infrared flux method. ICasagrande et al. 112010 ). deriving the mass of the secondary 



usually requires an observation of the reflex motion of the primary. Given the large number of EBs 
present in transit surveys, the ability to select a subset that have the precise masses of interest for 
the project is more efficient. 



- 4 - 



2. The Mass-Radius Relationship of K and M Dwarfs 



One example of binary science that benefits greatly from the merging of Doppler and photomet- 
ric survey data is the study of the stellar mass-radius relationship. Observations of EBs have been 
in excellent agreement with theoretical ni odels with the exception of K and M dwarfs (M < 0.8 Mq, 
Lopez-Morales 1120071 : iTorres et al. Il2010l ). Stars in this mass range have partially convective outer 
atmospheres, and can have radii that are 10-20% larger than models predict. Stars a bove this mass 



range have primarily radiative atmospheres and agree w ell with current models fe.g.. iPopper 1119971: 



Imbert II2002I : iTorres et al. II2008I : iMeibom et al. 



20091 . and additional references in ITorres et al. 



One possible cause of this inflation is the stronger magnetic f ields for K and M dw arfs in 
short-period binaries as compared to isolated stars of the same mass (IChabrier et al. 1120071 ). Tidal 
synchronization increases the rotational period of the stars to match the orbital period, and the 
stronger magnetic dynamo acts to suppress convection in the outer layers of the star. This results 
in an increase in cool starspots, reducing the effective temperature, and therefore the radius of 
the star expands to maintain the luminosity of the star. Recent theoretical work has shown that 
starspots predominantly distributed at polar latitudes, and/or suppression of convective transport 
in the outer layer s of the atmosphere, a re both effects caused by magnetic fields that can reproduce 



the observations (jMorales et al. Il20ld ) 



In order to place constraints on the interaction of magnetic fields and the stellar radius, a 
large sample of EBs with masses and radii measured to precisions of a few percent is required. 
Specifically, the sample should have at least one component in each EE that is 0.3 < M < 0.8 Mq 
and should have a wide range of orbital periods in order to sample a range of magnetic field 
strengths. However, despite several decades of surveys for photometric variability, there remains a 
paucity of such EB systems with precisel y determined masses a nd radii. For example, there were 
only 13 EB systems in the compilation of iLopez-Morales I (j2007l ) that fit the mass criterion above. 
Combining photometric and Doppler survey data can help increase the sample size by factors of 
several over the next few years, particularly at longer orbital periods (P > 5 days). Ground-based 
surveys for exoplanets identify many SBs and EBs that are, by design, bright stars (V < 12). This 
translates to higher signal-to-noise, requiring less telescope time for follow-up, and therefore a larger 
number of these binaries can be measured with high pr ecision. The Kepler rn i ssion is uncoverin g 



thousands of EBs with exquisite photometric precision (jCoughlin et al. I l2010l : iPrsa et al. I l2010l ) 



but a majority of those binaries are several magnitudes fainter than EBs found from ground-based 
exoplanet surveys, resulting in resource-intensive spectroscopic follow-up. In addition, any such 
ground-based observations are in direct competition with programs dedicated to exoplanet follow- 
up. 

Although several groups are working on studying the mass -radius relationship of double-lined 



spectroscopic binaries (SB2s) from ground-based surveys (e.g.. lDevor et al. II2008I). less attention 



is placed on binaries consisting of F/G primaries and K /M secondaries (but see 



Fernandez et al. 
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2009 ) ■ These binaries are generally brighter than K/M pairs, and avoid any biases introduced by 
studying binaries with mass ratios close to unity. The difficulty with studying binaries that have 
small flux ratios (-F2 is that the mass of the fainter component usually cannot be measured 

using the survey data. In order to measure the Doppler motion of the fainter companion, large- 
aperture telescopes equipped with high resolution sp ectrographs must be used in conjunction wit h 



two-dimensional correlation techniques (TODCOR. IZucker &: Mazeh Ill994l : IZucker et al. II2003I ). 
Only a minimum number of epochs per target using these large telescopes are required, because 
other parameters such as the orbital period, semiamplitude of the primary, and epoch of transit 
are known from the survey data. Observations can be taken near the RV extrema to measure the 
semiamplitude of the fainter component and thereby obtain both dynamical masses. 

As a case study, we present initial results from a cross-referencing of 24 SBs with P < 20 days 
from th e MARVELS Pil ot Project (hereafter MPP) with photometry from the Super WASP transit 
survey (jPollacco et al. This cross-referencing was performed after the RV observations 



from the MPP were completed, and the SuperWASP lightcurves were visually examined to identify 
which of the SBs were eclipsing. Out of those 24 SBs, five did not have a sufficient number of 
photometric observations and two are blended with other stars. Out of the remaining 17 SBs we 
have found two EBs: TYC 0272-00458-1 and TYC 1422-01328-1, which we refer to as TYC 272 
and TYC 1422 throughout the rest of the paper. TYC 272 is a single-lined spectroscopic binary 
(SBl), but we constrain the mass of the primary using several different techniques to estimate the 
properties of the host star. We then calculate the mass of the secondary, as well as the radii of 
the primary and secondary, from the survey data. TYC 1422 is found to be a triple system with 
all three components visible in high-resolution spectra. We note that although this work analyzed 
the photometry and spectroscopic data separately, many software pack ages exist to simultane ously 



model all observational data to derive the best-fit binary solutions, (e.g.. lPrsa &: Zwitter II2005I ). We 
demonstrate the success of cross-referencing Doppler and transit survey data and lay the foundations 
for more comprehensive studies using significantly larger data sets. 



3. Survey Data 



3.1. MPP Doppler Observations 



The MPP was a trial survey conducted in 2007 that used t he Keck ET instrum ent (jCe et al. 



2006al ) on the Sloan Digital Sky Survey (SDSS) 2.5m telescope (jOunn et al. 1120061 ) at the Apache 



Point Observatory. Although the primary objective was to serve as a test study for the MARVELS 
extrasolar planet survey, the MPP obtained between 5 and 38 RV measurements of 708 stars 
over 1-5 month baselines. The 59-object, fiber-fed instrument observed a total of twelve fields, 
and successfull y demonstrated the ability to conduct a multi-object, DFDI (dispersed fixed-delay 



2002l : lGe II2002I : lErskine II2002I : lErskine et al. Il2003l l survey for extrasolar 



interferometry, iGe et al 
planets. The spectral resolution of the Keck ET was R 



5100 and had a wavelength range of 
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495 < A < 585 nm. Despite its relatively small scale, the MPP was also able to demonstrate the 
ancillary science that can be done with such a survey, including the discovery of a sho rt-period, 
brown dwarf candidate with minimum mass rn-min = 64.3 Mjupiter (jFleming et al. 1 120101 ) . 



3.2. SuperWASP Photometry 



The WASP instruments provide flux measurements for millions of stars using wide-angle images 
of the night sky over a band pass of 400-700 nm defined by a broad-band filter. The survey 
has produced more than 35 transiting exoplanets to-date, o rbiting stars of A through K spectra l 



types and with orbital periods . 78 < P < 8.1 days (e.g., ICollier Cameron et al. 



Hebb et al. 


2010 


; ( 


Dueloz et al. 


2010) 


2007; 


Antipin et al. 


2010: 


Wils 


2010). 



2007a, 



201C; 



2010|), in addition to studies of variable stars (jNorton et al. 
2010l ). Eight cameras on each instrument provide images covering 
approximately 7.8x7.8 degrees using Canon 200mm f/1.8 camera lenses and e2v 2048x2048 CCDs. 
Synthetic aperture photometry usi ng an aperture radius of 49 arcseconds at the position of cataloged 
stars is performed on the images (jPollacco et al. 1120061 ). We extracted 1003 observations from the 
WASP archive for TYC 272 from 10 nights with good quality data covering either the primary or 
secondary eclipse. For TYC 1422, 649 data points from 12 nights of good quality data were used. 



4. Follow-up Observations 

4.1. Hereford Arizona Observatory Absolute Photometry 

We used the Hereford Arizona Observatory (HAO), a private facility in southern Arizona (ob- 
servatory code G95 in the lAU Minor Planet Center) to measure multi-band, absolute photometry 
of TYC 272 and TYC 1422. The facility has been used to measure transiting exoplanets with high 
photometric precisioij^, and has recently been used to observe EBs and variable star lightcurves 
with similar levels of precision. HAO employs an 11-inch Celestron Schmidt-Cassegrain (model 
CPC 1100) telescope fork- mounted on an equatorial wedge and located in a dome. The 11-inch 
telescope's CCD is an SBIG ST-8XE CCD with a KAF 1602E detector. HAO also includes a 
14-inch Meade Schmidt-Cassegrain (model LX200GPS) telescope, fork-mounted on an equatorial 
wedge located in a dome. That telescope's CCD is an SBIG ST-IOXME with a KAF-3200ME 
detector. A 10-position filter wheel accommodates SDSS and Johnson/Cousins filter sets. 

The stars were observed with Sloan g',r',i' filters with the 11-inch telescope. To derive the 
target star's magnitude in band x we use a generic photometry equation: 

= - 2.5 log (^] - Kz + S^C (1) 



• ^ttp : //brucegary . net/book_EOA/x . htm 



-7- 



where is the desired magnitude for observed band x, M^g is an instrumental zero-shift for band 
X, Fx is the measured intensity in the observed band, t is the integration time, K' is the zenith 
extinction for the observed band, z is the airmass, Sx is the star color sensitivity for the observed 
band and C is a star color defined using two bands (e.g., {B — V) or ( p' — r')). For the Johnsori - 
Kron -Cousins bands , standard stars are taken from the list published bv lLandolt fc: Uomoto I (j2007l ) 
and (Landol> 



Smith et al. 



Smith et al. 



20091 ) . For the SDSS bands, standard stars are taken from the list published by 
2OO2I) . B and V magnitudes are determined using the conversion equations given in 
(I2OO2I ) that convert g' , r' and i' magnitudes to -B, ^ and Rc, while Ic is determined 



using the i' fluxes. 



4.2. SMARTS Echelle Spectroscopy 

We obtained spectra from the SMARTS 1.5m telescope at CTIO using the R ~ 42000 Echelle 
spectrograph that covers 402 < A < 730 nm. These observations were used to resolve the com- 
ponents of TYC 1422 and get individual RV measurements for each star in the system. Seven 
epochs were observed using ten minute integrations between February and April 2010 with typical 
signal-to-noise ratios of ~ 20 per resolution element at A = 600 nm. 

The spectra were processed in the standard way for cross-dispersed Echelle spectra, using a 
pipeline written specifically for spectra taken with the Echelle spectrograph on the CTI0-1.5m. 
The routine processes the data using biases and quartz lamp observations taken at the beginning of 
the night, median combines three individual images while performing cosmic ray rejection, extracts 
the individual orders from the combined image and performs the wavelength solution on each order 
using a ThAr arc lamp taken either before or after each set of science exposures. 

A total of 22 orders are used each night to derive radial velocities via cross-correlation with 
a standard template. We use the bright GO star HD 108510 as the heliocentric RV standard 
star. Each spectral order is cross-correlated separately, then an iterative 3-0" clipping is performed 
prior to performing a weighted average to obtain a final RV measurement for each night. The 
components are identified each night via the peak and width of each feature in the CCF. The 
typical RV precisions ranged from 0.3 to 1.7 km s~^ for the various components. 



4.3. APO Echelle Spectroscopy 



We obtained ARCES (jWang et al. II2003I ) spectroscopy from the APO 3.5m telescope to check 
for additional spectra from unresolved companions. ARCES is an Echelle spectrograph with spectral 
resolution R ~31500 that covers the entire wavelength range from 320 < A < 1000 nm. on a single 
2048x2048 SITe CCD. Three exposures for a total integration time of 18 minutes were obtained 
for TYC 1422 on 2009 March 23, while two 20-minute exposures were taken of TYC 272 on 2010 
June 19. The spectra were reduced using an IRAF script. Spectra are corrected for bias and dark 
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subtraction, cosmic rays and bad pixels. Flatfielding is done using a combination of two different 
sets of quartz lamp exposures; one with a blue filter in place and another set without the filter. 
These two different sets of flats allow for optimal extraction of both the blue and red orders. A 
ThAr lamp is used for wavelength calibration. The co-added spectrum for TYC 1422 had a signal- 
to-noise of ~ 75 per resolution element at 600 nm, while TYC 272 had a signal-to-noise of ~ 105 
per resolution element. 



Case Studies 



5.1. TYC 272 



5.1.1. Orbital Parameters 



We first identified the binaries spectroscopically using the MPP RV data. Following lFleming et al 



(|2010l ). the RV uncertainties for each star are scaled by a "quality factor" that is a multiplicative 
scaling that modifies the formal RV uncertainties to better account for systematics defined as: 



Q 



RMS (X- (X)) 
MEDIAN (cjx) ' 



(2) 



where X represents the RV measurements and RMS is the root-mean-square residual. The quality 
factor is derived from the RMS of the RVs for the other 58 objects in the same field-of-view as the 
target, under the assumption that a majority of the objects per field should not have significant RV 
variability. The field containing TYC 272 has a Q = 3.44, suggesting the formal uncertainties are 
underestimated. Table [T] contains the times of observations, RVs, and both the formal and scaled 
RV uncertainties for TYC 272. A total of 14 epochs were obtained with an average, scaled RV 
uncertainty of 0.345 km s~^ for TYC 272. 



TYC 272 was fit using the RVSIM software program (jKane et al. II2007I ). Fig. [U shows the 
best-fit orbital solution to the MPP RV data. The x^/dof, after scaling by Q, is 2.62 with eight 
degrees of freedom. Because it is a short-period binary, there is additional RV jitter expected due 
to chromospheric act ivity, however, chromospheric effects are typically on the order of tens of m s~^ 
(jSantos et al. Il2000l ). The eccentricity is consistent wit h a circularized orbit, as expected for the 
best-fit orbital period P = 5.7282 days (jMazeh I bood ). The semi-amplitude K ~ 54.6 km s"^ 
corresponds to a minimum mass of the secondary mmin ~ 0.64 Mq assuming the mass of the 
primary is M = 1 Mq, making this EB an excellent system to further study in the context of the 
mass-radius relationship. The best-fit orbital properties are listed in Table [2j 
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5.1.2. Stellar Properties 



We utilize a spectral synthesis technique as part of the p rocess to estimate basic st ellar param- 
eters such as Tefj, log (5) and [Fe/H] for the primary. In IValenti &: Fischer I (j2005l ). they out- 
line a general procedure for deriving stellar properties using Spectroscopy Made Easy (SME; 



bod), 



Valeiiti &: Piskunov 1 1 19961 ). We follow the exact methodology described in IValenti &: Fischer 



but include some changes that allow for a more in-depth search of space. Specifi 



cally, we use the same spectral regions, line lists and continuum regions as in IValenti &: Fischer 



(j2005l ). but do not allow any of the individual abundances to be free parameters. 



Our SME spectral analysis pipeline is divided into two steps, because the spectral synthesis 
parameters are highly correlated, and therefore a range of initial guesses is required to adequately 
explore the space. First, a forward modeling procedure is used to calculate values for model 
spectra over a large grid of parameter space. This allows us to focus our analysis on the region of 
parameter space with the most likely, global, best-fit stellar parame ter solution. Second, we u se the 
SME Levenberg-Marquardt minimization algorithm as was used in IValenti &: Fischer I (120051 ) to fit 
observed spectra with synthetic spectra for 100 different initial input parameters. We select these 
different combinations of input values from the 100 grid points with the lowest chi-square values 
determined from the first step of our analysis pipeline. Accessing Vanderbilt Universitys ACCRE 
parallel computer facility allows us to efficiently explore the dependence of these initial parameters 
on our output parameters, as well as providing a good measure of the internal precision of SMEs 
minimization procedure. 

Using our implementation of SME on the ARCES spectrum of TYC 272, we derive best-fit 
values of Tog = 5905 ± 177 K, log (5) = 4.83 ± 0.19, metallicity [Fe/H] = -0.28 ± 0.12. The 
value for log {g) is larger than expected for stars of this temperature, but we note that log {g) is 
a very challenging parameter to measure spectroscopically, even at high spectral resolution, and is 
particularly challenging at the ARCES spectral resolution of ~ 31500. 

In addition to SME spectral synthesis, we analyzed the spectrum following a standard spectro- 
scopic method based on the requirements of excitation and ionization equil ibria. Our automate d 



analysis merges the approaches described in lPorto de Mello et al. I (j2008l) and 



and will be described in a forthcoming paper of the MARVELS series iFleming et al 



Ghezzi et 



mm 



2011 



). A 

total of 57 Fe I and 3 Fe II lines were used, after a 2-a clipping was applied to remove those lines 
with too high or too low iron abundances. The equivalent widths (EWs) were automatically mea- 
sured using the task "bplot in IRAF. We have derived a T^^ = 5720 ± 79K, log {g) = 4.20 ± 0.10 
and metallicity [Fe/H] = —0.37 it 0.10. We investigated the effects of a contaminating spectrum by 
modifying the measured EWs by a correction factor equal to the inverse of the expected flux contri- 
bution from the primary. In the case of a solar-type dwarf with a mid-K companion, the expected 
flux ratio is ~ 10% in V band, and the EW correction factor is ~ 1.1. The general effect is to 
increase the derived effective temperature, surface gravity and metallicity. However, the increase is 
small and well within the internal errors of our technique. We therefore conclude that a flux ratio 
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of 10% or lower would not significantly affect our measured values based on the EW method. 



We a,lso examine the Ha line profile , which serves as a T^s and chromospheric activity indicator, 
following iLvr a &: Porto de Mello I (|2005l ) and derive a best-fit T^s = 5461 it 118 K. This temperature 
is lower than the value determined via the Fe line EW method, but is affected by contamination 
from the secondary that fills in the Ha wings, resulting in an effective temperature determination 
that is systematically cooler. An analysis of the Hq core filling suggests that TYC 272 is much more 
chromospherically active than the Sun, however, there is a degeneracy between flux contamination 
from a companion and core emission due to activity, both of which can fill in the core of the Ha line. 
We further note that it is often the case that Ha-based and color-based Tcfj are systematically lower 
than temperatures determiri ed spectroscopically, one possible cause being chromospheric activity 
(jPorto de Mello et al. II2008I ). Given the relatively short orbital period of 5.728 days, it is likely 
that some amount of tidal spin-up has occurred, resulting in stronger magnetic fields and increased 
starspot activity. 

We complement the spectroscopic analysis with a variety of additional techniques to estimate 
stellar parameters based on other obser vables. To further estimate the luminosity class, we use a 
reduced proper motion (RPM) diagram (jCollier Cameron et al. Il2007bl ) to determine whether the 
host star is a likely giant or a dwarf /subgiant. The RPM diagram is excellent at distinguishing giants 
from dwarfs/subgiants. We tested the RPM technique on the Kepler Input Catalog (KICjfl, which 
has had its determ inations of luminosity class verified using Kepler's high-precision photometry 



(jKoch et al. Il20ld ). We find that only 2% of KIC giants were classified as dwarfs/subgiants using 



our RPM analysis. We use the proper motion values from the Tycho-2 catalog and NIR photometry 
from 2MASS (Table E]) to place TYC 272 on the RPM-J diagram (Fig. H bottom). The sohd, 
curved line in the RPM diagram is the dividing line between giants and dwarfs. The labels show 
the approximate regions where dwarfs, subgiants and giants are located. Given how precise the 
RPM diagram is at identifying giant stars, the fact that TYC 272 lies in the dwarf/subgiant region 
of the RPM diagram means it is highly unlikely that TYC 272 is a giant star. 



We use the infrared flux method (IRFM) (ICasagrande et al. Il2010l ) to estimate the effective 
temperature of the primary star via a Monte Carlo approach. We sample the observed V, J and 
Ks magnitudes assuming Gaussian errors, sample uniformly in [Fe/H] from —1.0 < [Fe/H] < 0.5 
and uniformly in extinction from < Ay < 0.0 75, where the upper bound is determined from the 
Schlegel et al. dust maps (jSchlegel et al.lll998l ). Fig. [2] shows the results of the Monte Carlo IR 
flux method (top). We estimate the effective temperature based on the (V-Ks) color to be T^ff = 
5459 (+233, —208) K, where the uncertainties are 1-a equivalents centered on the median of the 
asymmetrical T^g distribution. The difference in temperature between the IRFM and spectroscopic 
results can be caused by flux contribution from the secondary. Indeed, the peak of the ( J-i&)-based 
temperature distribution is slightly cooler, which can occur because the flux ratio in the NIR bands 
is larger than in the optical bands. This IRFM T^s is therefore a lower bound on the true T^s of 



^http : //eir chive . stsci . edu/kepler/kiclO/search.php 
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the primary due to reddening by the companion. 



To further study the effective temperature of the primary, we make use of the spectral energy 
distribution (SE P). We use the opti cal band measurements from HAO, near-IR measurements 



via the 2MASS (jSkrutskie et alJl2006l ) Point Source Catalog and UV mea surements from GALEX 



( Martin et al. Ibood ) . which are compiled in Table El NextGen models from Hauschildt et al. (jl999l ) 



are used to construct the theoretical SED. Based on the results from the spectroscopic and IRFM 
analysis, we construct a model with T^s = 5785, log (g') = 4.5 and [Fe/H] = —0.3 and compare 
with the observed fluxes. These stellar parameters were selected as representative values between 
the spectral synthesis, EW measurements. Ha profile analysis and IRFM temperature calculations. 
Distance and extinction are free parameters 

Fig. [3] (top) shows the result of the SED fit and the model spectrum. The blue points are the 
passband-integrated model fluxes for each filter, while the red crosses are the bandpasses of each 
filter (horizontal bars) and the measured error in the fluxes (vertical bars). A clear excess in the 
near-IR is seen when we do not consider flux from the secondary. Adopting a primary mass of 1 Mq 
and a radius of 1 Rq, we calculate a secondary mass of 0.64 Mq (K5 spectral type dwarf). We 
then add flux contribution from a second star with T^q = 4400K and a radius of 0.7 Rq, and find 
that the inclusion of flux from a K-type companion can explain the near-IR excess seen in the SED 
Fig. El (bottom). This SED analysis does not provide an independent measurement of the stellar 
parameters, but rather serves as a plausibility experiment that demonstrates flux contribution from 
a K-type companion can explain the apparent NIR excess if the primary has a T^g ~ 5785 K. 

Because the secondary is too faint to have its spectral features detected in the MPP spectra, 
the mass o f the primary canno t be directly measured via RV. Instead, we make use of Padova 



isochrones ([Marigo et al. II2008I ) to estimate the mass of the TYC 272 primary. Fig. U] plots radii 
and masses as a function of effective temperature for two assumed ages (100 Myr and ~9 Gyr) and 
a metalicity Z =0.01 ~ [M/H] = —0.28. This metallicity is chosen based on the SME and EW 
determinations, which are in good agreement with each other. The uncertainty in [M/H] does not 
significantly affect the derived range of probable masses for TYC 272, therefore we elect to use a 
single value for the tracks' metallicities, for clarity. The derived Teff values from the SME, EW, Ha 
and IRFM are plotted as horizontal bars representing the l-cr confidence intervals. We elect to use 
a Tcfj range equal to the EW 3-a confidence interval (5483 < T^s < 5957 K) as our best estimate 
for the primary Teff because it is expected to be minimally affected by flux contribution from the 
secondary. We note however that all four T^s are consistent to within 2-0", and in fact the SME, 
EW and IRFM temperatures all agree to within 1-a. Based on this T^fj range, we find the Padova 
models constrain Mi = 0.92 ± O.IMq for log (g) > 4.0 and 8.0 < log (age) < 9.95. 
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5.1.3. Light Curve Parameters 

The lightcurve for TYC 272 is shown in Fig. [5l with the best-fit model overplotted as the soHd 
hne. The top panel contains the primary eclipse, while the shallowness of the secondary eclipse 
can be seen in the bottom panel. The secondary eclipse of TYC 272 appears to be grazing, lead- 
ing to model degeneracies between the ratio of the radii and the surface brightness ratio. This 
is a well-known pro blem for grazing EBs, and the solution is to constrain the flux ratio when fit- 
ting the lightcurve (lAndersen et al. Ill983l). Th e flux ratio is usuall y constrained spectroscopically 



(|Clausen et al. 1 1200,4 ISouthworth et al. 1 120041 : Blake et al. 1120081 ^. and we utilized the ARCES 



spectrum to search for any evidence of a secondary spectrum. The ARCES spectrum was ob- 
tained at an orbital phase corresponding to an expected velocity shift of -27.25 km s^^, however, 
the signal-to-noise was not sufficient to detect a secondary spectrum for al l flux ratios. Fig. El 



compares the ARCES s pectrum with the (lower-resolution, R=12,000) FOE (jMontes et al. Ill99 



I 



Montes &: Martin Ill998l ) spectrum of 9 Cet, a G2 V standard. Two different wavelength regions are 
examined that contain absorption features present in stars of mid-F through mid-M spectral types. 
No set of secondary spectra are seen, and we estimate an upper limit of < 20% for the ffux ratio 
at these wavelengths. In the absence of sufficiently high resolution spectra, we estimate the ffux 
ratio using the mass range expected for Mi in §5.1.21 and the measured RV semiamplitude of the 
primary Ki. The ffux ratio is constrained to be 8 it 4%, the estimated ffux ratio in 1/ of a G-type 
primary and K-type companion. As shown in the SED analysis, such a ffux ratio reproduces the 
observed SED well, particularly in the 2MASS NIR bands. 

To model the lightcurves, we used the iktebofS implementation of the EBOP lightcurve model 



(|Nelson fc Davis I Il972l : IPopper fc Etzel I Il98ll ) to produce a least-squares fft to the data. We 



modified jktebop to include the normalization for each night of data as additional free parameters. 



We consult ed various tabulations of linear limb darkening coefficients (Ivan Hamme I Il993l : I Claret 
2000, l2004l ) and decided to adopt a limb darkening coefficient of 0.6 in the WASP band for the 



G-type primary of TYC 272. For the K-type secondary star we used a linear limb darkening 
coefficient of 0.7. An uncertainty of ±0.1 is assumed for these parameters, and the small effect 
of this uncertainty is accounted for in the standard errors quoted for other parameters of the 
lightcurve solution. The values and uncertainties of e cos uj and e sin oo from the RV orbital solution 
are included as constraints to the least-squares fit to the lightcurve. Without including such a 
constraint, the best-fit value of esinw = 0.2 it 0.07, however, because this is determined via the 
widths of the eclipses, we do not get a robust estimate from the lightcurve data alone. 

Table m summarizes the geometric parameters derived from fitting the photometry. To account 
for correlated errors ("red noise") in the photometry, uncertainties in the parameters are derived 
using a "prayer-bead" bootstrap Monte Carlo method. This method uses the residuals from the 
best fit to create synthetic lightcurves with noise characteristics that are similar to the real data 



' ^http : //www . astro . keele . ac .uk/~jkt/codes/ jktebop. html | 
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via circular permutation. Comparison of the parameters that are common to both the RV analysis 
and the photometric analysis show excellent agreement. The orbital period agrees to within 1-cj, 
while the ecosw and esinw values are consistent with a circular orbit. 



5.1.4- Determination of Mass and Radius 

The mass of the secondary, as well as the radii of the two binary components, can be determined 
in an SBl if a primary and secondary eclipse are both observed and the mass of the primary star 
is known (or assumed). This is accomplished via the mass function: 



(Ml +M2f GP^ 



(3) 



Where G is the gravitational constant, P is the orbital period, and oi is determined observa- 
tionally using the definition of the RV semiamplitude Ki from the radial velocity equation: 



«, = %^^J = (1.2) (4) 

2 vr sm [z) 

where e is the orbital eccentricity, i is the line-of-sight orbital inclination and the subscript j = {1,2} 
represents the primary and secondary, respectively. The radii of the two stars can be calculated 
directly from the geometric lightcurve parameters Ri and i?2 found in Table H] and the 
semimajor axis a = oi + 02. 

The uncertainties for Mi, M2, Ri and R2 are determined as follows. The uncertainty of Mi is 
taken from the independent measurement of the primary's mass, if such a determination exists, or 
from the range of estimates for Mi based on the star's color and/or spectra. Using implicit partial 
differentiation, the propagation of error equation, and removing the cross terms since Mi and C 
are independent, the uncertainty in M2 from Eqn. [3] is given by: 

_ Ual^C^ (Ml + M2)^ + al, (Mi + M2)n 
""^^ \ 3(M2)'-2C(Mi + M2) J 

where C is the right-hand-side of Eqn. [3) 

and ac is found by using the propagation of error equation and assuming P, Ki, e and i are all 
independent. Likewise, the uncertainties in the stellar radii are calculated via the propagation of 
error equation using the measured uncertainties in i?i a^^, R2 a^^ (Table Hj) and the calculated aa 
found via propagating the uncertainties in oi and 02. 

We adopt a primary mass of Mi = O.92ibO.lM0 based on the stellar characterization presented 
in gSTDand calculate M2 = 0.610 ±0.036 M©, Ri = 0.932 ±0.076 Rq and R2 = 0.559 ±0.102 Rq 
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(Fig. [7]). Both the primary and secondary are consistent (l-cx) with the model predictions for mass 
and radii. Metahcity is a minor effect on the mass and radius of the isochrones, so only the tracks 
corresponding to Z = 0.01 are shown. 

Although the relative errors for the masses and radii in the specific case of TYC 272 are 
larger than typically used to study the mass-radius relationship (~ 6% for M2, 8% and 18% for 
Ri and R2, respectively), relative errors of a few percent are possible if better constraints on 
the primary mass could be determined and more precise photometry of the grazing secondary 
eclipse co uld be obtained. Specifi c ally, large-aperture t elescopes and multi-dimensional correlation 



software (IZucker &: Mazeh I Il994l : IZucker et al. I l2003l ) can be used to extract RV solutions for 



faint companions, such that the masses of both components can be measured without relying on 
models or empirical relations to constrain Mi. In addition, high precision photometric follow-up 
of the shallow secondary eclipse will improve the precision in the radius. TYC 272 is not an 
ideal benchmark system because of the shallow, grazing secondary eclipse, however, the analysis 
presented here demonstrates the required steps to derive precise masses and radii for EBs using a 
combination of transit and RV exoplanet survey data. A summary of these analysis steps for EBs 
with low flux ratios, such as TYC 272, are described in ^ 



5.2. TYC 1422 

5.2.1. Orbital Parameters 

TYC 1422 had a total of 16 epochs observed during the MPP, with a, Q = 3.07 for its field 
and an average, scaled RV uncertainty of 0.304 km s~^. When fitting the MPP RV measurements 
for TYC 1422, the residuals were found to be exceedingly large (RMS ~ 6 km s~^). Such a large 
scatter is often caused by contamination of an additional spectrum from an unresolved source, a 
likely scenario for short-period SBs with similar fiuxes. 

The MPP DFDI pipeline was not designed to treat multiple sources within a single spectrum, 
and can produce large systematic errors in those cases. It is important to note that the inability 
to obtain precise RV measurements for SB2s applies to the MPP instrument only. Other survey 
instruments, including the MARVELS survey instrument itself, will be able to treat the case of 
double- lined binaries. The MARVELS survey pipeline stores information about the fringe fitting 
that can be used to remove the fringes from the stellar DFDI spectrum and re-create a "traditional" , 
de-fringed stellar spectrum. This spectrum can then be used to perform cross-correlation analysis 
and identify and measure multiple CCF components, thereby measuring RVs of double-lined binary 
stars. 

Because the MPP instrument passed starlight through an Iodine cell, a similar reconstruction 
of the stellar spectrum is not possible due to the Iodine lines that are present. Attempting to model 
the Iodine lines and remove them from the spectrum prior to de- fringing is also difficult, because 
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the Iodine lines are numerous and very narrow. In fact, most of the Iodine hues are unresolved 
at the resolution of the MPP instrument. For this reason, independent RV follow-up observations 
(SMARTS data) were obtained to measure the RVs of the multiple components. In general, such 
additional follow-up is not required and RV measurements for both components of double-lined 
binaries can be obtained directly from the survey data itself. 

Despite the large scatter in the residuals, the derived orbital period from the MPP data is 
in good agreement with the period derived from the Super WASP photometry (~ 6.2 days). To 
derive an improved orbital solution, we use the higher resolution SMARTS spectra to resolve the 
components and measure the Doppler reflex velocities of the individual components. We do not use 
the MPP data in our final analysis because it is so heavily affected by systematic errors. However, 
for completeness, we present the MPP RV data in Table O 

The TYC 1422 system consists of three sources, as demonstrated in the triple set of Na D 
doublets from the ARCES spectrum in Fig. [8j We do not attempt an SED analysis due to the 
complexity of the system, however, the basic stellar properties, including the HAO absolute pho- 
tometry, are compiled in Table [6l The RVs for each component are phase-folded to the orbital 
period and ephemeris as determined by the SuperWASP photometry in Fig. [9j The RV measure- 
ments for each component are connected by different line types for visualization. Components 2 
and 3 are identified as the eclipsing binary pair. Their RVs are in antiphase, as expected for an 
orbiting pair, and have relative velocities that cross the zero-line near the phases of primary and 
secondary eclipse. Component 1 is found to have a gradual linear trend of ~ — 75 m day~^. 
Fig. [lO] shows the RVs from Component 1 and the best-fit linear relation. The large uncertainty 
for the third epoch is caused by blending. Table [7] contains the times and measured RVs for each 
of the three components. 

The SMARTS RVs were fit using MPRVFIT2, an IDL procedure written by N. De Lee which 
combines a periodogram searching algorithm and a non-linear least squares fitter to find the optimal 
Keplerian orbit for a given set of radial velocit y points. The per iod search algorithm consists of 



a modified Lomb-Scargle period ogram ba s ed on ICumming I (120041 ) with an analytical False Alarm 



Probability (FAP) derived from Baluev (j2008l ). The software allows for simultaneous fitting of 
RVs from both components of a binary system. The SMARTS spectra consist of seven epochs, 
resulting in only one degree of freedom for a full Keplerian fit. However, because Components 2 
and 3 are a physically bound pair, we can exploit the fact that several parameters are identical 
for the two components. Specifically, the period (P), eccentricity (e), epoch of periastron (Tp) 
and systemic RV (7) are identical, while the argument of periastron (cj) differs by vr radians. This 
results in seven degrees of freedom between the 14 combined data points, the seven parameters being 
{P, e, uj, Tp, Ki, K2,j}, where Ki and K2 are the RV semi-amplitudes of the two components. Due 
to the linear trend observed in Component 1, we also include a linear term that reduces the degrees 
of freedom to six. Fig. [TT] shows the best-fit orbital solution for Components 2 and 3, and the 
orbital parameters are presented in Table El We find a best-fit period of 6.2005 days and a slightly 
eccentric orbit with e = 0.16 and a mass ratio of 0.78. The orbital period is in agreement with the 
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photometric period to within 1 — a. The slope of the hnear term is —48 it 18m s~ day" , which is 
of comparable magnitude and direction as the slope measured on Component 1, and therefore may 
be residual instrument drift. 



5.2.2. Light Curve Parameters 



We performed the same analysis using the jktebop software as done for TYC 272, however, 
third light (flux contributed from Component 1, L3 in Tabled]) is included as a free parameter in 
the least-squares fit to the lightcurve of TYC 1422. In principle one can compare a spectroscopic 
flux ratio based on the EWs of spectral lines with the best-fit flux ratio derived from the lightcurve 
solution. However, in the case of TYC 1422 and the ARCES spectrum, most of the lines are 
blended at the resolution of the spectrograph. The Na D lines shown in Fig. |8] are quite visible, 
however, the wings of the features are still blended, and it is challenging to find an appropriate 
pseudo-continuum from which to measure the EWs. Furthermore, the EWs of the Na D doublets 
are known to depend on T^q (jTripicchio et al. 1119971 : iDiaz et al. 1 120071 ) . and the Teg of the third 
component is not constrained in our analysis. 

We use limb darkening coefficients of 0.55 for the primary and secondary. The third light is 
accounted for within the WASP band, which ranges from 400-700 nm. and is similar to g' + r'. 
The lightcurve with the best-fit model overplotted is shown in Fig. [T2j Table [9] summarizes the 
geometric parameters derived from photometry fitting. Correlated errors are accounted for using 
the "prayer-bead" bootstrap Monte Carlo approach, as described in ^5.1.31 Timing of the primary 
and secondary eclipse events confirm the system is in an eccentric orbit. 



5.2.3. Determination of Mass and Radius 

Since both components of the EB pair have Doppler measurements, we can measure the masses 
and radii of Components 2 and 3 directly from the RV semi-amplitudes and the lightcurve. We find 
the semimajor axes oi and 02 using the RV semiamplitudes from Eqn. HI The masses can then be 
determined using the center of mass equation and the complete version of Kepler's Third Law: 

where Mi and M2 are the masses of the primary and secondary, respectively. 

The radii of the two stars can be calculated directly from the geometric parameters measured 
from the lightcurves in Tableland the value of the semimajor axis a = 01+02. The primary 
star is evolved off the main sequence and has a mass of Mi = 1.163 it 0.034 Mq and a radius of 
Ri = 2.063 ± 0.058 Rq. The secondary has a mass of M2 = 0.905 ± 0.067 Mq and a radius of R2 = 
0.887 lb 0.037 Rq and is consistent with a main-sequence G-type dwarf. The relative uncertainties 
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for Ri and R2 are 2.8% and 4.2%, very close to the desired precision goal of < 3%. Additional, high 
precision photometry of the eclipses will be helpful to further improve the precisions. The relative 
uncertainties for Mi and M2 are 2.9% and 7.4%. The larger uncertainty for M2 is in part due to 
the larger uncertainties in the RV measurements for that component. The SMARTS spectra had 
relatively low signal-to- noise (~ 20 per resolution element), and obtaining higher signal-to-noise 
spectra will reduce the uncertainty in M2 to better than 3%. Fig. [13] compares the masses and 
radii with Padova isochrones at an age of log (t) = 9.77 and metalicities of Z = {0.01, 0.019, 0.03}. 
TYC 1422 is consistent with solar metalicity (solid line), whereas the dashed lines are additional 
metalicity tracks for comparison purposes. 



6. Framework For EBs From Surveys 

When the RV semiamplitudes for both components can be measured spectroscopically (SB2s), 
the masses and radii for both stars can be determined directly from observable quantities in the 
survey data. TYC 1422 is such an example, although it is not a good benchmark system due to the 
presence of a third spectrum. Additional spectra can blend with the other components and affect 
their derived RVs, and the additional flux contribution complicates the photometry modeling. 
Despite these challenges, we were able to obtain 2.9% relative precision on Mi, 2.8% relative 
precision on Ri and 4.2% relative precision on R2. The MPP data were unable to be de- fringed 
due to the presence of the Iodine lines from the gas cell, and therefore follow-up observations 
with SMARTS data were required. Obtaining better signal-to-noise spectra would reduce the 
uncertainties in the masses further. We note that the MARVELS survey instrument does not pass 
starlight through an Iodine cell, and spectra of SBs with similar flux ratios can be de-fringed and 
analyzed using cross-correlation techniques like the ones used on the SMARTS spectra of TYC 
1422 without additional follow-up. 

EBs with smaller flux ratios require model-dependent constraints on the primary mass to derive 
the mass of the secondary and both radii, until observations with larger telescopes can be obtained. 
In this paper, we provide a framework for identifying EBs of interest from SBls using survey data 
and catalog information. We summarize the steps taken in a suggested "analysis flow diagram" 
presented in Fig. [TH This framework applies only to those EBs that have both a measurable 
primary and secondary eclipse. Once the mass of the primary and an uncertainty are adopted 
via comparison with isochrones and/or using empirical relations, the mass of the secondary, both 
stellar radii, and their uncertainties can be calculated using the equations presented in §5.1.41 

We apply this technique to TYC 272, however, the grazing secondary eclipse makes it a poor 
benchmark EB as well. The 5.9% relative error on M2 is dominated by the large uncertainty for the 
estimate of Mi . Observations with high res olution spectrographs a t high signal-to- noise can be used 



to derive Ki and K2 using TODCOR (see iBender &: Simon 



20081 . for an example in the near-IR). 



Such a technique can exploit the fact that the other orbital parameters are well-determined from 
the survey data, and only a measurement of K2 is required to directly measure Mi and thereby 
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complete the analysis. In addition, obtaining a high signal-to-noise, high resolution spectrum at a 
phase where the spectral lines are well-separated can yield spectroscopically-determined flux ratios, 
which are then used as additional constraints in the lightcurve modeling of eclipses that are grazing. 

In favorable cases where there are deep, total eclipses and no third li ght contamination, trans it 



survey photometry itself can yield radii to better than 2% relative error (jSouthworth et al. II2OIII ). 
However, in both SB2 and SBl situations, the survey photometry may sometimes be insufficient 
to obtain such high precisions. Because the epoch of central transit is well-determined from the 
survey data, follow-up observations can be planned when a primary or secondary eclipse is known 
to occur. In addition, obtaining transits in multiple filters can aid in the EB analysis. We note 
that this can be a problem for any study of EBs using transiting exoplanet survey data, and is not 
unique to the methodology presented here. 



7. Conclusion 

We have presented a framework for eclipsing binary science via cross-referencing of photometric 
and Doppler exoplanet survey data. We have found two new EBs around F/G primaries; TYC 
272, an SBl consisting of a K dwarf secondary, and TYC 1422, a triple-lined SB that includes 
an evolved primary star and a G dwarf secondary. Although neither of these two specific binaries 
are ideal benchmark systems for studying the mass-radius relationship, we have shown that the 
combination of survey data can be a rich source of EBs, for which many physical parameters can 
be directly measured or constrained from the survey data itself. In particular, the analysis of TYC 
272 provides a framework for estimating the masses and radii of SBls from survey data. Obtaining 
additional, high resolution, high signal-to-noise spectra can then directly yield the semiamplitudes, 
and thus masses, of both components. Such an approach can be used in particular to study the 
mass-radius relationship of F/G + K/M EBs. 
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±cLDic 1. 


MPP RV Observations - 


TVP 979 
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RV 


(Tuv (formal) 
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(km s-l) 
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2454101.98701 
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0.086 
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2454195.81803 
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0.119 
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2454215.66891 


-16.141 


0.086 


0.297 


2454221.71415 


-6.716 
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0.339 


2454224.70505 


69.756 
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0.364 


2454249.66233 


-19.993 


0.100 


0.343 


2454251.67471 


65.708 


0.101 


0.349 


2454255.67428 


-18.045 


0.145 
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Table 2. Orbital Properties - TYC 272 



Parameter Value l-cr Uncertainty 



Period (days) 


5.7282 


0.0003 


K (fems-l) 


54.584 


0.294 


e 


0.001 


-HO. 009 
-0.001 


Ul 


5.548 


0.006 


Tp 


2454102.9539 


0.0074 


Sys. vel. 70 (km s~^) 


34.430 


0.231 




2.62 





Table 3. Stellar Properties - TYC 272 A 



Parameter Value l-tr Uncertainty 



fia (mas yr 
W (i 



a (J2000)'' 


176.54255245 (deg) 


11:46:10.21 


(HH:MM:SS) 


S (J2000)'' 


1.68753726 (deg) 


+01:41:15.13 


(DD:MM:SS) 




21.443 




0.113 


NUV'^ 


15.398 




0.0042 


B 


10.945 




0.035 


V 


10.266 




0.018 


Rc 


9.849 




0.022 


Ic 


9.460 




0.035 


vl 


12.055 




0.100 


a' 


10.547 




0.010 


r' 
i' 


10.045 
9.880 




0.025 
0.020 


z' 


9.800 




0.026 




8.893 




0.021 




8.520 




0.036 


-f'''52MASs'^ 


8.443 




0.023 


COS ((5))^ 
mas yr~-^)^ 


146.9 




1.6 


-138.3 




1.6 


RPMj 


5.417 







='Tycho-2 Catalog l|H0g et al. Ibood 'l 
t^GALEX jMartin et al. Il2005f) 

'^2MASS l lSkrutskie et al.ll2006l) Point Source Catalog 



Table 4. Light Curve Properties - TYC 272 



Parameter 


Value 


l-cr Uncertainty 


J (Surface brightness ratio at disk center) 

-R1+-R2 

a 


0.14 
0.096 


0.01 
0.004 


R2 

Ri 

i (degrees) 


0.60 
86.6 


0.12 
0.3 


e cos uj 
e sin a; 
Period (days) 
To (HJD) 


0.0004 
-0.002 
5.72840 
2454534.7020 


0.0009 
0.001 
0.00001 
0.0007 


«i 

a 


0.060 


0.002 


R2 

a 


0.036 


0.006 


^ (400-700 nm band) 
e 

UJ (degrees) 


0.048 
0.002 
284 


0.026 
0.0008 
101 
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Table 5. MPP RV Observations - TYC 1422' 



JrlJJJuTC 


-tlV 


cr {formal) 


cr jiY (scaled) 




(km s~"^) 


(km s~'^) 


(km s""*") 


24^41(11 Sfi'^Tfi 


10.848 


0.090 


0.276 


2454102.03097 


-12.295 


0.128 


0.393 


2454105.97331 


67.537 


0.109 


0.335 


2454128.81685 


49.394 


0.099 


0.304 


2454136.77639 


62.745 


0.092 


0.282 


2454136.81223 


53.464 


0.093 


0.286 


2454164.75574 


-11.161 


0.106 


0.325 


2454165.75527 


33.872 


0.090 


0.276 


2454165.79120 


33.716 
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0.279 


2454186.65694 


52.251 
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0.279 


2454191.72814 


53.452 


0.099 
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2454194.74076 


13.140 


0.090 


0.276 


2454195.73280 


-11.254 


0.090 
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2454221.62566 


34.133 


0.094 


0.289 


2454224.61948 


33.477 


0.094 


0.289 


2454254.63340 


59.181 


0.129 


0.396 



^These data are dominated by flux contamination and are not 
used in the final RV analysis. 



Table 6. Stellar Properties - TYC 1422 System 



Parameter 


Value 




l-CT Uncertainty 


a (J2000)'' 


152.41204974 


(deg) 


10:09:38.89 (HH:MM:SS) 


S (J2000)'' 


17.58113689 


(deg) 


-fl7:34:52.09 (DD;MM:SS) 


9' 
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0.010 


r' 
i' 




9.797 
9.726 


0.008 
0.013 


^2MASs'' 




8.875 


0.026 


-'^2MASs'' 




8.641 


0.028 


^^•52MASs'' 




8.563 


0.023 


fici (mas yr~^ cos (S))^ 
tJ-S (mas yr-l)" 




-30.6 


0.9 




-5.6 


1.0 



''Tyclio-2 Catalog l lH0g et al. Ibood) 

''2MASS jSkrutskie et al.ll2006l) Point Source Catalog 
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Table 7. SMARTS RV Observations - TYC 1422 





Component 1 


Component 2 


Component 3 


HJDuTC 


RV 




RV 


<^RV 


RV 






(km s-l) 


(km s-l) 


(km s-i) 


(km s-i) 


(km s-l) 


(km s-i) 


2455236.67830 


43.107 


0.298 


-18.420 


0.602 


132.754 


1.650 


2455241.64797 


41.238 


0.450 


67.401 


0.601 


9.110 


2.507 


2455254.64277 


31.593 


10.625 


31.593 


10.625 


72.102 


1.688 


2455270.64841 


40.024 


0.488 


90.128 


0.490 


-12.069 


1.648 


2455278.59449 


38.709 


0.901 


78.638 


0.655 


0.463 


1.490 


2455284.60109 


39.142 


0.304 


87.278 


0.607 


-8.814 


1.602 


2455293.53225 


38.567 


0.381 


-10.746 


1.071 


108.512 


2.509 




Table ! 


3. Orbital Properties - TYC 1422 








Component 2 




Component 3 


Parameter 


Value 


1-a Uncertainty 


Value 


l-cr Uncertainty 


Period (days) 


6.2005 




0.002 


6.2005 


0.002 




K (km s-l) 


65.366 




0.516 


84.046 


1.297 




ecosw 


-0.1494 




0.0123 


0.1494 
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esincj 


0.0576 




0.0170 


-0.0576 


0.0170 




Tp 


2455236.8133 




0.0554 2455236.8133 


0.0554 


Sys. vel. 70 (km s ^) 


45.026 




0.356 


45.026 


0.356 


Linear Trend d (km 


s ^ day ^) 


-0.048 




0.018 


-0.048 


0.018 




xVdof 


4.135 






4.135 





Table 9. Light Curve Properties - TYC 1422 



Parameter 


Value 


l-o" Uncertainty 


J (Surface brightness ratio at disk center) 


0.88 


0.01 


R1+R2 
a 


0.162 


0.002 


R2 

Ri 


0.43 


0.03 


i (degrees) 


87.7 


0.5 


ecosoj 


-0.1283 


0.0003 


esinw 


0.092 


0.011 


L3 (Contamination from 3'^'' Comp.) 


0.25 


0.08 


Period (days) 


6.199450 


0.000007 


To (HJD) 


2454194.3663 


0.0008 


Hi 
a 


0.114 


0.003 


«2 

a 


0.049 


0.002 


^ (400-700 nm band) 


0.16 


0.02 


c 


0.158 


0.006 


oj (degrees) 


144 


3 
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TYC 0272-00458-1 




0.0 0,2 0.4 0.6 0.8 1.0 

Orbitol Phose 



Fig. 1. — Best-fit orbital solution for TYC 272 from the MPP data. The RV uncertainties have 
been scaled by a factor of 3.44 to account for additional systematic errors. 
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^ TYC 0272-00458-1 



0.2 0.4 0.6 0.8 1.0 

(J-H) 

Fig. 2. — Top: Monte Carlo calculation of effective temperature using the infrared flux method. We 
estimate a T^s = 5459l;208^ from this analysis. Bottom: RPM-J reduced proper motion diagram. 
The dividing line between giants and dwarfs is the solid curve. Approximate regions where dwarfs, 
subgiants and giants inhabit are labeled. The location of TYC 272 is consistent with a dwarf or 
subgiant. 



Av=0.075 T=5785 logg=4.5 [Fe/H]=-0.3 
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Fig. 3. — SED fits using NextGen models to TYC 272. Top: Model that only includes flux from the 
primary, assuming stellar parameters from the combined spectroscopic analysis. Bottom: Model 
that includes flux contribution from a K5 dwarf secondary, demonstrating the near-IR excess can 
be explained by a K-type companion. 
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Fig. 4. — Radius (top) and mass (bottom) as a function of effective temperature for different 

Padova isochroncs and metalicity Z = 0.01. The solid line is for an assumed young age (100 Myr) 
while the dashed line is for an assumed old age {^9 Gyr). The horizontal, solid lines are the T^s 
1-0" confidence intervals for the various stellar parameter techniques that were applied. Only stars 
with log (g) > 4.0 are shown, since there are no indications the primary has a lower log (g). We use 
these isochrones to constrain the mass of the primary (Mi) in order to derive M2, R\ and i?2- 
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Fig. 5. — Phase-folded SuperWASP lightcurve and best-fit model for the primary (top) and sec- 
ondary (bottom) eclipses of TYC 272. 
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9 Cet (top, G2.5V), TYC 0272-00458-1 (bottom) 




0.41 , , , , I 

5850 5855 5860 5865 5870 5875 



Wovelength (Angstroms) 

Fig. 6. — Comparison between the R ~ 12,000 FOE spectrum of the G2.5V standard 9 Cet (top 
of each plot) with the R ~ 30,000 ARCES spectrum of TYC 272 (bottom of each plot). The two 
wavelength ranges contain absorption lines present in stars mid-F though mid-M. No secondary set 
of spectra can be seen. 
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Fig. 7. — Radii and masses for the two components of TYC 272. Both stars are consistent (l-cr) 
with model predictions. Only a metalicity of .Z^ = 0.019 is shown, since it has only a minor effect 
on the masses and radii of dwarfs and subgiants. 
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TYC 1422-01328-1 




5886 5888 5890 5892 5894 5896 5898 5900 

Wovelength (Angstroms) 

Fig. 8.— ARCES spectrum of TYC 1422 centered on the Na D doublets. The three sets of doublets 
can be clearly seen, and components are marked with the same names assigned to them while 
measuring their RVs with the SMARTS data. The components were identified in the SMARTS 
CCF's by examining the width and shape of their CCF peaks. 
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TYC 1422-01328-1 - RVs and LC 




0.0 0,2 0.4 0.6 0.8 1.0 

Orbltol Phose 

Fig. 9. — Photometry and RVs for the three components of TYC 1422 phase-folded to the best-fit 
period and epoch of central transit in Table [9l The connecting lines between the measured RVs 
for each component serve as a visual aid. Components 2 and 3 are identified as the EB pair. 
Their measured RVs are antiphased as one would expect for an orbiting pair. They also phase 
appropriately such that the RV zero-crossing occurs near primary and secondary eclipses. Note the 
connecting lines do not cross at the precise primary and secondary eclipse phases because they do 
not represent the actual Keplerian orbit. 
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Fig. 10. — RVs for Component 1 in TYC 1422. The large error bar for the third epoch is caused 
by blending. The data are well-fit with a gradual slope of ~ —75 m s""^ day~^, which could be 
residual instrument drift given that a similar slope is found when fitting Components 2 and 3. 
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Fig. 11.— Best-fit orbital solution for Components 2 and 3 of TYC 1422 from the SMARTS 
observations. The black data points and curve are the RVs for Component 2 and the best-fit 
orbital solution, while the grey points are for the less massive Component 3. The jagged nature 
of the best-fit curve is a result of phase-folding the small linear-trend included in the fitting. The 
bottom panel shows the residuals to the fit. 
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Fig. 12. — Phase-folded SuperWASP lightcurve and best-fit model for the primary (top) and 
secondary (bottom) eclipses of TYC 1422. 
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Podovo Isochrones, Z = {0.01, 0.019, 0.03} 




Moss (solor mosses) 



Fig. 13. — Mass-Radius diagram for TYC 1422, Components 2 and 3. Padova isochrones are 
overplotted using a log (t) = 9.77. This age was selected as the best-fit by visual inspection of 
tracks with differing ages. We find excellent agreement using a solar-metalicity of Z = 0.019. For 
comparison, we also show tracks with Z = 0.01 (left track, dotted line) and Z = 0.03 (right track, 
dotted line). 
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Fig. 14. — Suggested analysis flow diagram that demonstrates the steps one would take to obtain 
masses and radii for an EB with a small flux ratio, starting with exoplanet survey data and catalog 
information. A dashed line represents a "no" decision. The analysis depicted here only works for 
EBs that undergo both primary and secondary eclipses. 



